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We report experimental measurements of the flow in a cubic container submitted to a longitudinal 
libration, i.e. a rotation modulated in time. Velocity fields in a vertical and a horizontal plane are 
measured in the librating frame using a corotating particle image velocimetry system. When the 
libration frequency 00 is smaller than twice the mean rotation rate flo, inertial waves can propagate 
| in the interior of the fluid. At arbitrary excitation frequencies Co < 2f2o, the oscillating flow shows 

(— i i two contributions: (i) a basic flow induced by the libration motion, and (ii) inertial wave beams 

propagating obliquely upward and downward from the horizontal edges of the cube. In addition to 
these two contributions, inertial modes may also be excited at some specific resonant frequencies. We 
characterize in particular the resonance of the mode of lowest order compatible with the symmetries 
of the forcing, noted [2,1,+]. By comparing the measured flow fields to the expected inviscid inertial 
modes computed numerically [L.R.M. Maas, Fluid Dyn. Res. 33, 373 (2003)], we show that only a 
(—| 'I subset of inertial modes, matching the symmetries of the forcing, can be excited by the libration. 

i 

I. INTRODUCTION 

Rotating fluids support the existence of a singular class of waves called inertial waves,— ~— which arc 
r/} ' anisotropic and propagate because of the restoring nature of the Coriolis force. These waves exist only for 

excitation frequencies <jq lower than twice the rotation rate In a confined volume of fluid, inertial waves 
' may become phase-coherent and experience a resonance due to multiple reflections over the container walls, 

leading to the so-called inertial modes. These inertial modes are relevant to geophysical and astrophysical 
flows (in liquid cores of planets and stars), but also to mechanical engineering flows (e.g., in liquid- filled 
projectiles). 

Inviscid inertial modes are the eigenmodes of a given container geometry, and can be found in general 
when the walls are either normal or parallel to the rotation axis^^ but also in some specific cases such 
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as spheres, spheroids^ and, to some extent, in spherical shells (namely for so-called R-modes that lack 
radial displacement). However, in general, when sloping walls are present, the wave focusing and defocusing 
induced by the peculiar reflection law of inertial waves^ precludes the existence of inviscid eigenmodes, and 
the concentration of energy along particular beams leads to wave attractors^ The resonance frequencies of 
inviscid inertial modes can be derived analytically only in some specific geometries, such as cylinders, spheres 
and spheroids liajfi In the case of a parallclcpipedic container, such as used in the experiments presented here, 
the frequencies and the spatial structure of the inviscid inertial modes have been determined numerically^ 
In geophysical and astrophysical flows, several types of forcing may be at the origin of the excitation 
of inertial modes. Inertial modes have been excited experimentally in a spherical cavity by a longitudinal 
. libratioQafl — i.e. a time modulation of the rotation rate — , and have recently been described in numeri- 

cal simulations in spheres and spherical shells^ - — Precession^— and periodic deformation of the walls 
modeling gravitational tide o 18 ' 19 are other examples of forcing, providing an efficient generation of inertial 
modes. These inertial modes have been proposed to contribute, through nonlinear self-interaction, to the 
generation of steady zonal flows, which are visible for instance in the atmosphere of gaseous planets like 
Jupiter *2r— 

Although axisymmetric geometries have been primarily considered to observe and characterize inertial 
modes, non-axisymmetric geometries such as parallelepipeds 9 - are also of fundamental interest. In partic- 
ular, the symmetries of the container, the symmetries of the forcing, and the role of viscosity, are critical 
parameters in determining which inertial modes can be excited. 
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I Visualization 

T 




2L = 30 cm 

FIG. 1. (Color online) Schematic view of the (30 cm) 3 water tank mounted on the rotating platform. The bottom 
and top walls are located at z = and z — 2L, and the origin of the coordinates (a;, y) is taken at the center of 
the square cross section. Two-dimensional PIV measurements are achieved in the vertical plane y = and in the 
horizontal plane z = 4L/3 in the librated frame using a co-rotating laser sheet and a camera aiming normally at it. 

Incrtial modes may be present in any laboratory experiment performed in a rotating container. They 
have been for instance detected in ensemble averages of turbulence generated by the translation of a grid in 
a rotating container^ - — When present, these modes may couple to the turbulence, and have a profound 
influence on its statistical properties. In particular, turbulence in the presence of inertial modes, which are 
non homogeneous by nature, cannot be considered as freely decaying, raising the issue of the relevance of 
the homogeneous framework to describe rotating turbulence in laboratory experiments. 

In this paper, we investigate the spatial structure of the oscillating flow generated by the longitudinal 
libration of a cube, and we characterize the efficiency of this configuration to excite incrtial modes. Mea- 
surements are performed both in a vertical and in a horizontal plane (the rotation axis is vertical) using a 
co-rotating two-dimensional particle image velocimetry (PIV) system. For libration frequencies lower than 
twice the rotation rate oq < 2Qo, we observe, in addition to the basic oscillating flow induced by the libra- 
tion, the propagation of oblique inertial wave beams emitted from the horizontal edges of the cube. These 
wave beams, similar to the ones observed in cylindrical geometry^2£ originate from the convergence of 
Ekman fluxes near the edges of the container. In addition, for a specific libration frequency in the explored 
range ao/2£lo e [0.60; 0.73], libration also excites an inertial mode, of spatial structure in close agreement 
with the inviscid computation of Ref. [^. Our results suggest that, among all the cigen modes predicted 
numerically, only a small subset of low order modes matching exactly the symmetries of the cube libration 
can be actually excited. 

II. EXPERIMENTAL SETUP 

A. Flow geometry and rotating platform 

The experimental setup, sketched in Fig. |TJ consists in a closed cubic glass tank, of inner size 2L — 30 cm, 
filled with water and mounted on a precision rotating turntable of 2 m in diameter. We use a coordinate 
system centered on the horizontal square, —L < x,y < L, with the horizontal walls at z = and z = 2L, 
perpendicular to the rotation vector f2 = f2(t)e z . The mean angular velocity of the rotating platform is 
set to fio = 0.419 rad s _1 (4 rpm). The Ekman number of the system, which compares the viscous to the 
Coriolis force, is E = v/(2VIqL 2 ) = 5.3 x 10 -5 . The rotation of the fluid is set long before the start of 
the libration, at least 1 h, in order for transient spin-up recirculations to be damped. Once the solid-body 
rotation is reached, the turntable is submitted to a longitudinal libration, which consists in a modulation of 



3 



the angular velocity Q(t) around f2o a t a frequency <jq with a peak-to-peak amplitude 2ef2o> 

n(t) = n [1 + e cos(ct *)] . (1) 
In the frame rotating at constant velocity fig , the libration of the cube is described by the angle 

(p(t) — e— sin(crot). (2) 

The normalized libration amplitude e, which is the Rossby number of the problem, is varied between 2% 
and 16%. In order to excite inertial waves, the libration frequency cto can take in principle any value in the 
range [0,2il ]. In this paper, we focus on the restricted range <t /2^o G [0.60; 0.73]. For these frequencies, 
the peak-to-peak libration angle Aip = 2ef2o/co lies in the range from 1° to 13° for e = 2% to 16%. The 
relative precision in the control of the instantaneous rotation rate is better than 10~ 3 . After the start 
of the libration, we wait at least half an hour before the data acquisition (which represents 8te, where 
te = L^vQq)^/ 2 is the Ekman timcscalc) in order to reach a stationary regime. 



B. Particle image velocimetry (PIV) measurements 



Velocity fields are measured in the librated reference frame using a two-dimensional PIV system mounted 
on the rotating platform (Fig. [T]) . 26 ' 27 Measurements are performed either in the vertical plane y — or in 
the horizontal plane z = 4L/3 = 20 cm. This last particular choice is motivated by the fact that z = 4L/3 
does not correspond to a node for the inertial modes considered in this paper (the centered horizontal plane 
z = L is a node for the inertial modes of even vertical wavenumber n). The flow is seeded with 10 /xm tracer 
particles, and illuminated by a corotating laser sheet generated by a 140 mJ Nd:YAG pulsed laser. For both 
horizontal and vertical measurements, the entire 30 x 30 cm 2 flow sections are imaged through transparent 
sides of the tank with a high resolution 2048 x 2048 pixels camera aiming normally at the laser sheet. 

Each acquisition consists in at least 1000, and up to 3000, images taken at a sampling rate between 12 co 
and 36 oo which correspond to at least 80 libration periods. The sampling rate is chosen according to the 
libration amplitude in order to keep a typical particle displacement of the order of 5 pixels between two 
successive images. PIV fields are computed over successive images using 32 x 32 pixels interrogation windows 
with 50% overlap, leading to a spatial resolution of 3.5 mm. This resolution is not enough to resolve the 
thickness of the Ekman boundary layers, Se = LE 1 / 2 ~ 1 mm, but is appropriate for the flow structures 
associated to inertial waves and modes in the bulk of the fluid. 



III. BASIC LIBRATION FLOW 



A. Inviscid solution 



We first characterize the basic flow response to the libration forcing. Since in the present experiments 
the modulation period 2tt/o-q ~ [20; 25] s is much shorter than the Ekman time te — 230 s (which is the 
timescale required for the angular velocity of the fluid to match that of the boundaries) , the core of the flow 
can be considered as essentially inertial, and rotating at the constant rotation rate . In the frame rotating 
at Qq, the flow is therefore approximately at rest, and surrounded by oscillating recirculations induced by 
the periodic motion of the walls. 

The inviscid response to the libration of an arbitrary container, whose walls are either parallel or normal 
to the rotation axis, can be derived by assuming that the flow is strictly two-dimensional. In the absence of 
viscosity, the absolute vorticity of the fluid in the frame of the laboratory, 

ui a = u> + 2Sl(t)e z , 

must be conserved, and hence given by the mean vorticity 2J7oe z . Here u> is the relative vorticity, as 
measured in the libration frame. If viscosity is present, this result remains approximately valid far from the 
boundaries, with the additional assumption that ao/Oo 3> E 1 ^ 2 (i.e., for a rapid libration compared to the 
Ekman timescale). With the total angular velocity of the libration given by Eq. ([T]), the relative vorticity us 
is therefore vertical, homogeneous, and given by 

w,(t) = -jr 1 - - V" = - 2en o cos(a t). (3) 
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(a) Inviscid theory (b) Experiment 




FIG. 2. (Color online) Velocity fields excited by libration in a container of square cross-section, (a) Inviscid solution, 
computed from Eq. (|A3|I . (b) PIV measurements at the phase Got = (maximum libration velocity). The libration 
frequency is cro/2fio = 0.648 (corresponding to the mode [2,2,+], which is not resonant in this experiment), and its 
amplitude is e — 0.04. 

In an axisymmetric container, the resulting flow is an oscillating solid body rotation of angular velocity 
— eilo cos(<7oi): this simply describes a fluid at rest in the frame rotating at constant rate fV The case of 
a non-axisymmetric container is more complex, and can be solved in terms of a streamfunction rp(x, y, t) = 
^(x, y) cos(<7oi), such that (u x ,u y ) = (—d y ijj,d x i(j). Equation ^ therefore takes the form of a Poisson 
equation for the spatial part of the streamfunction, 

A* = -2efl - 

The solution of this equation for a square domain, subject to the condition that the streamfunction vanishes 
at the boundary, is given in the Appendix, in Eq. ()A3|) . The streamlines are circular near the center of 
the container, as in a solid-body rotation, but they become gradually square near the boundaries (see the 
velocity field in Fig. [2^). 

B. Experimental measurements of the libration flow 

Figure [5] compares the measured libration flow obtained by PIV and the inviscid solution (|A3|) . The 
velocity field shown here is taken at the phase aot = corresponding to the maximum amplitude of the 
libration, when the flow in the librated frame is anticyclonic (i.e., uj z = — 2ef2o)- The libration frequency (To 
is chosen here far from any resonant frequency of inertial mode. The agreement between the measured field 
and the inviscid solution is excellent in most of the flow section, except near the boundaries. It must be 
noted that the vanishing of the velocity near the boundaries may be affected by the PIV resolution (about 
3 mm), which is of the same order as the expected boundary layer thickness. 

An interesting feature of Fig. &b) is the wavy shape of the experimental iso- velocity lines: this is a first 
indication that, in addition to the basic libration flow, the flow also contains a wave component. This 
additional component is further described in the next Section. 

For all the libration amplitudes investigated here (e = 2% to 16%), the fluid response to the forcing 
is found to remain essentially linear. This is demonstrated in Fig. [3J where the temporal energy spectra 
E{o) are shown for different forcing amplitude e. Here, the spectrum is defined as E(a) = (\u a \ 2 ) xy , where 
u a (x,y) is the temporal Fourier transform computed at each location, and () xy is the spatial average over 
the horizontal plane. The Fourier transform is computed over at least 80 periods of libration. 

All spectra show a narrow peak at the forcing frequency ctq, well above the white noise level of the PIV 
measurements. We note that, since we are only interested here in forcing frequencies ctq > higher 
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FIG. 3. Temporal energy spectra E(a) for a forcing frequency <7o/2fio = 0.675 (corresponding to the mode [2, 1, +]), 
for three libration amplitudes e = 0.04, 0.08 and 0.16. 

harmonics na^, if present, are beyond the upper limit 2f2 , and are therefore not governed by the dynamics 
of inertial waves. 

We can note that, for all forcing amplitudes e, the peak at <jq contains at least 97% of the total flow energy. 
The remaining energy is associated to secondary peaks at a — 0, a = arid higher harmonics (a > 2Slo) • 
The secondary peak at a = f2o, mostly visible at low e, corresponds to a residual fluid motion synchronized 
with the platform rotation rate. The peak at a = can be associated to the generation of a weak mean 
flow. This mean flow is probably a nonlinear response to the libration forcing, and is not investigated in the 
present paper. 

In order to improve the signal-to-noise ratio, we perform in the following a band-pass filtering of the 
velocity fields at the forcing frequency <jq. This procedure consists in filtering the Fourier transform u CT at 
all a except in a narrow region centered on ao of width 8a/2Vl( i = 0.024 (which contains the energy of the 
peak to within 1%), and computing the inverse Fourier transform of the resulting filtered field. 

IV. IIWISCID INERTIAL MODES AND SYMMETRIES OF THE LIBRATION FORCING 
A. Inviscid inertial modes 

The eigen frequencies and the spatial structures of the inertial modes in a parallelepiped whose sides are 
parallel or perpendicular to the rotation axis have been determined numerically for an inviscid fluid in Ref. [^. 
These predictions have been achieved from the numerical resolution of the eigcnproblcm defined from the 
inviscid linearized equations in a rotating frame. Any frequency <7q in the range [0, 2f2rj] corresponds to an 
inertial mode. Accordingly, an infinity of inviscid modes may be excited if the system is forced at a given 
frequency with a finite bandwidth. 

Following the notation introduced in Ref. l24l, we label the inertial modes as [n, m, s]. The first index n is 
the normalized vertical wavenumbcr, such that the velocity field u shows n recirculation cells in the vertical 
direction. With the container walls taken at z — and z = H, the velocity field has the form 



with ui = u x e x + u y e y (we restrict to H = 2L in this paper). Inertial modes are therefore stationary in the 
vertical direction, but their horizontal structure u(x, y, t) may be either stationary (the so-called "sloshing" 
modes) or propagating. The second index m enumerates the eigen frequencies of modes of sign s from the 
largest one (m = 1) down. Larger values of m essentially correspond to finer structures in the horizontal 
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plane. Since each mode m is expressed in terms of an infinite amount of horizontal Fourier modes, the 
index m is not directly related to a number of nodes as for the vertical index n. Finally, the sign s refers 
to the invariance of the mode with respect to rotation of angle 7r about the z axis. More precisely, a mode 
u(x, y, z, t) has symmetry s if 



The resonance frequencies of inertial modes [n, m, s] are increasing functions of n and, at fixed n and s, 
decreasing functions of m. This behavior essentially originates from the dispersion relation for plane inertial 

waves, 



where 8 is the angle between the wavevector k and the rotation axis, cos 9 = fc z /|k|. Low frequencies a are 
therefore associated to nearly horizontal k, i.e. to small k z and hence to small n, or/and to large k XiV and 
hence to large m. 



B. Symmetry properties and boundary conditions for a viscous fluid 

In the case of a viscous fluid, the no-slip boundary condition at the walls is expected to affect the spatial 
structure of the inertial modes. The spectrum and spatial structure of viscous inertial modes in a rotating 
parallelepiped have not been computed yet. It is therefore of first interest to check to what extent the inviscid 
modes found numerically could be recovered in a viscous fluid at finite Ekman number. Since viscosity damps 
preferentially the high order modes (i.e. the modes with large indices n and/or m), we expect the libration 
to force more efficiently modes of low order. This is straightforward for the vertical wavenumbcr n, which 
is naturally associated to a viscous damping proportional to vn 2 . No equivalent simple law exists for the 
horizontal index m, but we can similarly expect that lower m, associated to larger scales in the horizontal 
plane, will be favored in the presence of viscosity. 

From the symmetries of the boundary conditions, it is possible to anticipate which inertial modes are 
compatible with the libration forcing. In the frame rotating at constant velocity Qq, the angular oscillation of 
the top and bottom walls is described by the velocity u(x, y, z = 0) = u(x, y,z = 2L) = eflo cos(crot) (—ye x + 
xe y ) which, according to Eq. ([5]), has symmetry s = +. Moreover, the no-slip boundary conditions impose 
equal horizontal velocity at the top and bottom walls, which is satisfied only for even vertical wavenumbers 
n (the basic libration flow described in the previous section is vertically invariant, and is hence associated 
to n = 0). Finally, there is the possibility that particular symmetries associated to the second index m may 
even further reduce the set of modes compatible with the symmetries of the libration forcing. According 
to these symmetry properties, the longitudinal libration is expected to excite only a subset of the modes 
[n,m, s] among those having even n and s = +. 

It must be noted that those symmetry arguments state which modes are forbidden by this forcing, but 
they do not state which modes, among the allowed ones, are effectively excited by the libration. In other 
words, we would like to address the question of how the fluid motion induced by the oscillating walls, and 
in particular close to the edges of the cube, will be transmitted to the inertial modes, and to characterize 
the efficiency of this transmission. 

In the following, we investigate the flow response for libration frequencies in the vicinity of two inertial 
modes of low order allowed by the symmetries of the forcing: [2,1,+] and [2,2,+]. More specifically, we 
systematically characterize the flow in the range <jq/2£Io <G [0.60; 0.73] surrounding these two modes. The 
list of all modes in this range, truncated at n < 4 and m < 6, is given in Tab. Q] Among these modes, only 
the four ones [2,1,+], [2,2,+], [4,5,+] and [4,6,+] (marked by a * symbol) have the correct symmetries (even 
n and s = +) and could be observed in principle. However, we will show in the following that viscous effects 
and additional symmetry properties actually select only the single mode [2,1,+] in this list. 





a = 2Q,q cos 6, 
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mode 


frequency 


[n, m, s] 




[2, 1,+]* 


0.6742 


[2,2,+]* 


0.6484 


[3. 2, -] 


0.7271 


TO Q 1 

[3, 3, -] 


0.680/ 


To /i 1 

[3,4,-] 


0.6848 


f3 3 +1 
|o, -rj 


0.6485 


[3,4,+] 


0.6258 


[4,5,+]* 


0.6960 


[4, 6,+]* 


0.6945 


[4, 5, -] 


0.6889 


[4, 6, -] 


0.6643 



TABLE I. List of all the inviscid inertial modes [n, m, s], truncated at n < 4 and m < 6, in the range of frequencies 
<to/2S1o 6 [0.60; 0.73]. n is the normalized vertical wavenumber, m characterizes the horizontal structure, and s = ± 
refers to the symmetry or antisymmetry of the mode with respect to rotation of angle n about the rotation axis (see 
Ref. for details). The * symbol marks the 4 modes having symmetries compatible with the libration forcing: even 
n and s = +. 
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FIG. 4. (Color online) Spatial structure, in the vertical plane y = 0, of the oscillating flow excited with e = 0.02 and 
at frequencies <to/2Qq — 0.60 (a) and 0.73 (b), in the absence of inertial modes. The velocity fields are extracted 
by temporal band-pass filtering around the excitation frequency. The arrow fields illustrate the (x, z) components 
of the velocity field at a given arbitrary phase of the oscillation, and the gray shade/color maps the corresponding 
normalized vorticity uj y /max\ujy\. Resolution of the velocity fields has been reduced by a factor 5 for better visibility. 
In (a), bold arrows show the direction of the phase velocity for each wave beam. In (a) and (b), the solid black 
lines show the expected direction predicted for each wave beam from the dispersion relation, cos 6 = <ro/2f2o. 



V. EDGE BEAMS AND INERTIAL MODES 
A. Flow in the vertical plane: non-resonant case 

We first describe the flow in the vertical plane y = in the absence of inertial modes. The velocity u x ,z 
and the vorticity uj v (x, z) shown in Fig.@]for e = 0.02 are for forcing frequencies cto/2^o = 0.60 (a) and 0.73 
(b), chosen far from any low order inertial mode. Fields are not displayed in lateral strips of 12 mm width 
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along the walls due to poor quality of the data induced by light reflections. Since the basic libration flow is 
normal to the plane y = 0, the vertical flows here show the x and z deviations from the libration flow, and 
contain therefore only the contributions from the wave component of the flow. 

We observe that the libration excites a cross shaped pattern made of four oblique shear layers originating 
from the bottom and top edges of the cube. These shear layers actually correspond to propagative inertial 
waves in two-dimensional plane wave packets, which contain approximately one wavelength in their trans- 
verse direction. There are actually eight two-dimensional wave packets emitted from the four upper and the 
four lower edges of the container; only the four ones emitted from the edges defined by x = ±L, z = 0, 2L 
can be seen in the vertical cut y = 0. The direction of the beams is found in good agreement with the 
dispersion relation (J7J, as shown by the black lines in Fig. @] making an angle 6 = cos _1 (o"o/2fio) with the 
horizontal. In these wave beams, fluid particles describe anticyclonic circular motion at frequency <jq in the 
planes inclined at angles ±0 with the horizontal (the sign depends on the considered beam). The shearing 
motion traced by the vorticity uj y is related to the variation of the wave phase in the direction normal to ±9. 
We also observe that the phase of the wave travels normal to the beam, with a phase velocity c v directed 
towards the vertical wall from which originates the beam. In Fig. EJa) [9 = cos _1 (0.60) ~ 53°], the four 
wave beams show interferences at their intersections, which make the wave pattern rather complex. On the 
contrary, in Fig. [4](b) [9 = cos _1 (0.73) ~ 43°], the opposite wave beams are almost aligned for this particular 
frequency. As a consequence, the wave beams show constructive interferences, leading to almost standing 
waves along the diagonals, but propagative outwards on their sides. 

These edge beams arc similar to the ones observed in cylindrical containers, e.g. in the early experiment 
of McEwani (forced by a tilted top lid) and in the numerical simulations of Duguct et al^ (forced by 
an oscillating axial compression). Similar beams are also found in the recent simulations of Sauret et 
aL— (forced by a longitudinal libration), although in this case the beam angle is apparently not directly 
related to the libration frequency, and probably results from the turbulence in the boundary layers over the 
sidewalls. In all these references, the container is cylindrical, resulting in conical edge beams, emitted from 
the two circular edges and focusing towards the rotation axis. In our non-axisymmetric geometry, since the 
disturbance sources are linear segments, there are eight two-dimensional edge beams emitted from the 4 
upper and the 4 lower edges of the container. 

These edge beams originate from the oscillating outward (rcsp. inward) Ekman layers on the upper 
and lower walls associated with the prograde (resp. retrograde) part of the libration motion. Far from 
the lateral walls, where the streamlines of the libration flow are approximately axisymmetric, this flow is 
approximately radial, of order ef^oV x 2 + y 2 . In the classical axisymmetric spin- up problem, the horizontal 
flow in the unsteady (but slowly varying) Ekman layer is vertically deflected into Stewartson layers along the 
lateral wallsi However, in the rapidly oscillating situation examined here, the matching condition between 
the Ekman and Stewartson layers cannot be satisfied, and the flow in the Ekman layer detaches when it 
reaches the lateral walls, generating an oscillating shear layer in the bulk^ of characteristic velocity efloL. 
As a consequence, although the velocity scale of the edge beams is inertial, viscosity is responsible for this 
flow feature. 

The width of these edge beams can be compared to predictions from a viscous theory. Because of viscous 
damping, it is known that a wave beam spreads as I 1 / 3 , where I is the distance from the sourcei 29 ' 30 More 
precisely, the thickening of the full-width half maximum 8(1) of a self-similar plane beam emitted by a linear 
disturbance of small size can be computed^ 

with E = v I (2floL 2 ). According to this model, a wave beam emitted from an edge reaches the central area of 
the container after a distance £* = L/sva.6 (for 9 > 45°), where its width is 5(1*) ~ 6.84(sin0)~ 2 / 3 L E 1 / 3 ~ 
4.2 — 5.0 cm for the range of angles considered here. This prediction provides correct agreement with the 
observations in Fig. [U showing beams of thickness of order of 3 — 4 cm near the center. Note that the case 
9 ~ 45° is specific: the overlap of the upper and lower wave beams leads to a different scaling,— 8 ~ E 1 / 4 . 

B. Flow in the vertical plane: resonant case 

We still examine here the structure of the flow in the vertical plane, but now at frequencies where inertial 
modes are expected. We select the two frequencies <7o/2fio = 0.674 and 0.648, corresponding to the inviscid 
predictions for the modes [2,1,+] and [2,2,+] respectively (see Tab. HJ. The velocity and vorticity fields are 
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FIG. 5. (Color online) Spatial structure, in the vertical plane y = 0, of the oscillating flow excited with e = 0.04 
and at frequency <jo/2fio = 0.674 (a-c) and cro/20o = 0.648 (d-f). (a,b,d,e) correspond to the experimental fields 
extracted by temporal band-pass filtering around the excitation frequency. (c,f) show the corresponding numerical 
inviscid modes [2,1,+] and [2,2,+] expected at the same frequency. The libration phase is aot = for (a,d), and 
trot = 7r/2 for (b,e). The gray shade/color maps the normalized vorticit y fi eld u] y /max\uj y \. Resolution of the velocity 
fields has been reduced by a factor 5 for better visibility (see also Ref . [33) . 



shown in Fig.O at the phases aot = (a,d) and tt/2 (b,e) of the libration. In the frame rotating at constant 
velocity fio, the phase aot = corresponds to ip — and Q,(t) = ef2o (i.e., maximum rotation rate), whereas 
aot = tt/2 corresponds to ip maximum and f2(t) = [see Eqs. (CO)-©]. 

Looking first at the velocity fields (b) and (e), both frequencies show two oscillating cells in the vertical 
direction, as expected for a mode n = 2 (see also the movies available at Ref. HH) . These oscillating cells 
are mainly visible at the phase a$t = tt/2, when the angle ip is maximum and the libration motion vanishes. 
At this phase, the velocity fields clearly satisfy the s = + symmetry, with u x (— x, z) = —u x (x,z) and 
u z {— x, z) = u z {x,z). On the other hand, at the phase a^t = corresponding to the maximum libration 
velocity (a,d), the x and z components of the velocity associated to this mode n = 2 vanish and one recovers 
the propagative edge beams already evidenced in Fig. [4] 

Looking now at the vorticity field reveals finer structures of the flow. At phase aot = (a,d), the typical 
cross shaped pattern composed of four edge beams is found for both forcing frequencies. On the other hand, 
at phase tt/2 (b,e), the shape of the vorticity field is found to depend now on the forcing frequency. Whereas 
it keeps approximately its cross-shape structure for cto/20q = 0.648 (e), it shows four nearly circular extrema 
of alternate signs for <7o/2ilo = 0.674 (b). This indicates that, at aot = tt/2, the four recirculation cells 
dominate the vorticity for cro/2^o = 0.674 (b), whereas their vorticity is partially hidden by the strong edge 
beams for a /2il = 0.648 (e). 

In order to understand the different behaviors between the two frequencies, we compare now the experi- 
mental fields with the numerical predictions of the inviscid modes [2,1,+] and [2,2,+], shown in Figs. 03c) 
and (f), respectively (movies of the numerical modes are also provided at Ref. |32| ). Note that the phase 
of the numerical fields is arbitrary since they are not related to any specific forcing: the numerical fields 
in Fig. [5] have been simply chosen here at phases which provide a good visual correspondence with the 
experimental fields. Although the vorticity patterns of modes [2,1,+] and [2,2,+] look similar, with four 
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FIG. 6. (Color online) Spatial structure, in the horizontal plane z = 4L/3, of the wave component of the flow u' 
excited with e = 0.04 and at frequency <ro/2fio = 0.674 (a-c) and <7o/2f2o = 0.648 (d-f). Same layout as in Fig. [5] 



extrema of alternate signs, their velocity fields arc very different: the vertical velocity is maximum on the 
vertical axis for [2,1,+], but it is zero for [2,2,+]. The smooth pattern of ui y found experimentally in Fig.[5jb) 
matches actually well the prediction for [2,1,+]. On the other hand, the comparison fails for Fig.[5te), both 
for the velocity and vorticity fields, suggesting that the mode [2,2,+] cannot be excited by the longitudinal 
libration. 

We can interpret these observations as follows. For <jq/2SIq = 0.674, the flow is a supcrimposition of 
the pattern of propagative edge beams and a resonant stationary [2,1,+] mode. Whenever the libration 
angle is zero (i.e. at phases a^t = and 7r), the instantaneous amplitude of the [2,1,+] mode vanishes, 
and only the edge beams remain visible. Turning now to the case ao/2f2o = 0.648, the mode [2,2,+] is 
apparently not excited, and the edge beam pattern remains visible all the time. Interestingly, it seems that 
a flow reminiscent of the mode [2,1,+] is excited instead (see Fig. Efe)), although with a weak amplitude, 
probably because of the closeness of the eigen frequencies (less than 4%) of the two modes. It will indeed 
be confirmed in Sec. I VII that the resonance of the [2,1,+] mode spreads over a significant frequency range 
because of viscosity. 



C. Flow in the horizontal plane 

We finally turn to visualizations in the horizontal plane z = AL/?>, still for the two frequencies ao/2QQ = 
0.674 and 0.648. In order to focus on the wave component of the flow, we remove in the following the 
libration component, 

u'(x, y, t) = u(x, y, t) - ui ib (x, y) cos(<7 *), (9) 

where u(x,y,t) is the total velocity measured by PIV, and uub(x, y) = V x (^e z ) is the inviscid libration 
flow computed from Eq. (|A3[) . In Fig. El we show the resulting horizontal velocity u' x y and the corresponding 
vertical vorticity ui' z {x,y) fields for both frequencies, for the same phases <?ot = and n/2 as previously. 
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Here again, measurements are not displayed in lateral strips of width 6 mm along the walls because of the 
poor quality of the PIV data. 

For both frequencies, the horizontal flow shows a complex spatial structure which verifies the relation 
u' xy {—x,—y) = —u' xy {x,y), confirming the s = + symmetry of the flow. Both frequencies show a clear 
square pattern of vorticity extrema, which is the trace of the interferences between the 8 propagative 
inertial wave beams emitted from the 8 horizontal edges of the container. Note that since the libration flow 
has constant vorticity far from the boundaries, one has uj' z = lu z + 2e£lo cos(erot), so the vorticity patterns 
computed from the total flow and from its wave component are identical. The amplitude of the wave 
component in the bulk of the flow is typically a factor 10 below the libration velocity scale c^IqL, indicating 
a weak efficiency of the excitation of inertial waves by libration. 

In Figs.jnjc) and (f), we show for comparison the numerical predictions of the inviscid modes [2,1,+] and 
[2,2,+] in the horizontal plane. The mode [2,1,+] (c) consists in a nearly axisymmetric flow, in which each 
velocity vector describes an elliptic oscillation in the anticyclonic direction. During one period of the mode, 
the flow goes through the following sequence: cyclonic, centrifugal, anticyclonic, centripetal (see Ref. HU). 
The mode [2,2,+] (f) has a very different structure: it has four vortices located along the sidewalls, with 
a hyperbolic point in the middle. This pattern is essentially rotating as a whole around the rotation axis, 
with only weak deformations near the walls due to the non-axisymmetry of the container. 

The comparison between the numerical inviscid modes and the experimental measurements is revealing. 
The structure of the mode [2,1,+] can be easily recognized in the experimental data, superimposed to the 
square pattern of the interfering edge waves. For instance, the numerical mode chosen at the phase shown 
in Fig. |6tc) matches well the experimental flow at the libration phase a^t = [Fig. EJa)]. At this phase, 
the flow is normal to the plane y = 0, in agreement with the vanishing of the mode in the vertical plane 
observed at the same phase in Fig. [SJa). On the other hand, the structure of the inviscid mode [2,2,+] 
cannot be recognized in the experimental data at any phased confirming that this mode cannot be excited 
by libration. For both frequencies, the vorticity field is essentially dominated by the interfering edge waves, 
and is always different from the vorticity of the predicted inviscid modes [Fig.[5Jc) and (f)]. 

The absence of mode [2, 2, +] in the experimental fields may be understood from symmetry arguments. 
In addition to the s = + symmetry, which corresponds to an invariance under a rotation of n about the z 
axis, the [2,1,+] mode turns out to be also invariant under a rotation of 7r/2 about the z axis [Fig. H2c)], 
which exactly matches the symmetry of the libration forcing for a cube. On the contrary, the mode [2,2,+] 
has only the 7r symmetry, and it is antisymmetric with respect to the rotation of 7r/2 about z. Taking into 
account this additional symmetry explains why only [2,1,+] is excited by the libration and not [2,2,+], even 
when the forcing frequency exactly matches the predicted frequency. Excitation of [2,2,+] would require 
two adjacent walls to move in opposition of phase, which is not possible with a rigid container. 

These observations confirm that the libration forcing is able to excite only a subset of the possible inertial 
modes. Excitable modes [n, m, s] have even n and s = +, and must in addition be symmetric under rotation 
of 7r/2 about z, which is satisfied for [2,1,+] but not for [2,2,+]. 

VI. RESONANCE CURVE 

We finally describe the energy of the oscillating flow as the libration amplitude e and the forcing frequency 
<7q/2D,q are varied. First, we have varied the libration amplitude e for the two frequencies cto/2^o = 0.648 
and 0.674 corresponding to the (excitable) mode [2,1,+] and the (non excitable) mode [2,2,+]. The energy 
of the wave component u'(x, y, t) of the flow [see Eq. ©] is plotted in Fig. [7] as a function of e for these two 
frequencies. This energy is computed over the horizontal components of the velocity in the horizontal plane 
z = 4L/3, 

k h = (u' x 2 + u' y 2 ) xy , (10) 

where the overbar stands for time average. 

For both frequencies, the scaling of the energy is compatible with e 2 , confirming that the wave component 
of the flow essentially responds linearly to the excitation amplitude e in the range 0.02 — 0.16. The energy 
for the frequency corresponding to [2,1,+] is a factor 2.5 above that for [2,2,+]. Since we know that the 
mode [2,2,+] is not excited, the curve for [2,2,+] essentially corresponds to the interfering propagating edge 
beams. The additional amount of energy in [2,1,+] therefore corresponds to the resonant mode [2,1,+]. 
These curves suggest that, in the horizontal plane, the energy of the mode [2,1,+] and that of the edge 
waves are of the same order. 
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FIG. 7. Normalized kinetic energy kh/{0,aL) of the wave component of the flow in the horizontal plane z — 4L/3 
as a function of the excitation amplitude e, for libration frequencies corresponding to the [2,1,+] (o) and [2,2,+] (o) 
modes. The dashed lines show the best fits in e 2 . 




C7o/2f2o 

FIG. 8. (Color online) Normalized kinetic energy k v /(eQoL) 2 in the vertical plane y = as a function of the 
excitation frequency <7o/2f2o for an excitation amplitude e = 0.02. The vertical lines indicate the eigen frequencies 
of the inviscid modes of low order (only modes with n < 4 and m < 6 are shown). Solid lines correspond to s = — 
modes and dashed lines to s — + modes. The size of the vertical segments is taken proportional to 1/n. 



In order to characterize in more detail the resonance of the mode [2,1,+], we have also computed the 
energy in the vertical plane y = 0, 

k v = (ul+ul) xz . (11) 

The libration flow being normal to the measurement plane, it does not contribute to k v , so this definition 
essentially describes the energy of the wave component of the flow. Focusing on a single value of the 
excitation amplitude, e = 0.02, we have performed a systematic scan over the excitation frequency ero/2f2o 
in the range [0.60; 0.73]. Since the energy basically scales as e 2 , the energy k v in Fig. [5] is normalized now 
by (efloL) 2 , which represents the order of magnitude of the energy of the basic libration flow. 

The resonance curve shows a well-defined peak of energy k v ~ 3 x 10~ 2 (etloL) 2 , centered on er /2S!o = 
0.674 ± 0.005, and of width at half maximum Aer/2f2o ~ 0.02. The peak frequency is very close to the 
numerical prediction for the inviscid mode [2,1,+], <7q/2Qq = 0.6742, confirming that this mode is the only 
one which is significantly excited in the frequency range investigated here. We can note a slight shift of about 
<7o/20o ~ +0.002 between the predicted and the observed peak, which probably originates from a viscous 
detuning effect. Although this detuning has not been computed for the inertial modes of a parallelepiped, it 
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is essentially similar to the one found for the inertial modes of a sphered The frequency shift for the sphere 
has been actually determined analytically, and is of order of o~q/2VLq ~ (0.02 — 0.05)\/E for low order modes. 
This would give <7o/2f2o — 0.002 — 0.004 for E = 5.3 x 10 -5 , which is consistent with the shift observed here 
for the cube. 

The resonance peak is surrounded by a nearly constant plateau at k v ~ (5 ± 1) x 10 -3 (efio£) 2 , which 
accounts for the wave beams emitted from the edges of the container. This plateau confirms that the edge 
beams are always present in the flow, and are hidden only at the particular frequencies at which a mode 
is excited. It is worth noting that the width of the resonance peak is significantly larger than the width of 
the peak at <7o in the temporal energy spectra of the flow (see Fig. [3]) . This suggests that the quality of the 
resonance of the mode [2,1,+] is not fixed by the precision of the libration forcing, but is rather governed 
by the viscosity. This is confirmed by the fact that the resonance width Acr/2f2o ~ 0.02 is of the order of 
(f2o T £;) _1 — 0.01, where te = L/^Mq is the Ekman timescale. 

As expected, inertial modes which do not have the symmetries of the libration forcing (e.g. [3,4,+] or 
[3,3,+], with odd n), are not found in the resonance curve. More interestingly, the even and symmetric (with 
respect to rotation of it) mode [2,2,+] is not found either, confirming that the antisymmetry (with respect 
to rotation of 7r/2) of this particular mode is not compatible with the symmetry of the libration. However, 
it appears that at this frequency the system still lies in the far tail of the [2,1,+] resonance, which probably 
explains the residual flow structure associated to the mode [2,1,+] found in Fig.[3{c). We can conclude that, 
because of viscous effects and symmetry properties, the libration forcing selects only one single resonant 
inertial mode in the range [0.60; 0.73]. The slight asymmetry in the tails of the resonance curve might be 
due to the presence of the modes [4,5,+] and [4,6,+] at frequencies slightly larger than the [2,1,+] peak. 
It turns out that [4,5,+] has the ir/2 symmetry, whereas [4,6,+] has the tt symmetry only; it is therefore 
conceivable that a slight contribution of the mode [4,5,+] is responsible for the asymmetric shape of the 
resonance curve. 



VII. EXCITATION OF INERTIAL MODES BY THE EDGE BEAMS 

The present results raise the question of the relation between the inertial modes and the edge beams 
originating from the convergence of Ekman fluxes near the edges of the container. Interestingly, for forcing 
frequencies close to the frequency of low-order symmetric inertial modes [n, 1,+] with even n, the edge 
beams form a periodic ray pattern in the vertical plane (x, z). This is illustrated in Fig. H^a) and (c) for 
libration frequencies corresponding to the modes [2,1,+] and [4,1,+], showing one and two X-shaped patterns 
respectively. This figure suggests that periodicity of the rays leads to the formation of a standing wave, since 
energy can propagate along these rays either way. It motivates the use of a two-dimensional simplification, 
in which the planar wave beams emitted from the edges of the cube propagate in the vertical plane along 
rays, at an angle given by the dispersion relation ([7]). This 2D approximation does not apply near the 
corners, where the flow must be three-dimensional. But away from the walls, the edge beams essentially 
protrude in the bulk of the fluid, as if the container were infinitely long in the along-edge direction. 

We consider in the following the 2D ray orbits in a vertical plane originating from the edges of the cube. 
These orbits can be characterized by integers, denoting the number of reflections from vertical (i) 

and horizontal (j) boundaries. Here we count edge reflections as reflections on both the vertical as well as 
the horizontal wall. The orbits in Fig. G^b) and (d) therefore classify as (1,1) and (2,1) periodic orbits 
respectively. The angle which the beam makes with the horizontal is 8ij = tan _1 (j/i), yielding a normalized 
frequency 



similar to the expression for the eigenfrcquencies of the transverse modes of an infinite channel^ The 
frequencies are dense in [O,2f2o]- This spectrum is degenerate: patterns and have the same 

frequency if i'/i — j'/j (such degeneracy is not obvious for the inertial modes in the parallelepiped). The 
resulting ray pattern may be associated with a set of cells, with hyperbolic points at the intersection between 
rays. Restricting to ray patterns compatible with symmetries of the libration, this construction yields 2i x 2j 
cells [see dashed lines in Fig. |9^b) and (d)], corresponding to an n = 2i inertial mode structure. No similar 
connection can be made between index m and j from this simple geometrical construction, because m is not 
related to the number of cells in the horizontal direction. 
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FIG. 9. (a,c) Flow in the vertical plane y = corresponding to the modes [2, 1, +] and [4, 1, +]; (b,d) corresponding 
ray-tracing patterns (1,1) and (2,1). Figures are taken at the libration phase aot = 0, for which the mode has zero 
amplitude, so only the wave beams are visible. 



If we compare the frequencies of the periodic patterns (i, j) and the frequencies of the inertial modes 
[n,m, s], we obtain a reasonable agreement for n = 2i, m = 1, and s = + (sec Tab. [TTJ) . The discrepancy 
is 4.7% for [2, 1, +], and decreases for larger n. We can also note that, for j > 1, it is always possible to 
find, among the modes having n = 2i and s = +, a mode m > 1 such that er^- ~ <y nms to within a few 
%. For instance, for (i,j) = (1,2), one has <7i2/2f2o = 0.447, which is close to the frequencies of modes 
[2, 4, +] and [2, 5, +] (0.465 and 0.434, respectively). One explanation for the apparent ability of strictly 2D 
periodic orbits to predict eigenfrequencies of the 3D problem quite accurately may lie in the fact that while 
the eigenspectrum is dense, there are indications that the density of states (the amount of eigenfrequencies 
per frequency increment) may be nonuniform. Computations in long channels^ show that eigenfrequencies 
cluster preferentially around the eigenfrequencies of strictly transverse modes of an infinitely long channel, 
and there are indications that this also applies to the higher vertical modes (large n) in a cube. 

The coincidence between the frequencies of inertial modes and periodic beam patterns of similar spatial 
structure suggests the following mechanism for the excitation of modes in a libratcd container. A given 
libration frequency induces a set of inertial wave beams originating from the detachment of the oscillating 
Ekman layers at the edges. If the libration frequency is close to an inertial mode, the resulting edge beam 
pattern is nearly periodic, and has the same symmetries as the mode — this is at least the case for the two 
modes [2,1,+] and [4,1,+] shown in Fig. [5] This suggests that the edge beams can supply energy to the modes 
directly in the bulk of the flow, with the correct vertical periodicity. The libration flow, being associated to 
n = (vertical invariance), is unable to excite these modes other than viscously, through the edge beams, 
because it lacks the vertical structure of the inertial modes. This may be different in the presence of sloping 
boundaries, which enforce vertical motions, and may excite inertial waves inviscidly. We may conclude that 
inertial modes are excited with an inertial velocity scale eSloL, but this excitation requires the presence of 
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pattern 


Oij 


<Jij /2Qq 


mode [n, m, s] 




Difference 


(1,1) 


45° 


0.707 


[2,1,+] 


0.674 


4.7% 


(2,1) 


26.6° 


0.894 


[4,1,+] 


0.875 


2.1% 


(3,1) 


18.4° 


0.949 


[6,1,+] 


0.938 


1.2% 



TABLE II. Comparison between the frequencies of the periodic edge beam patterns (i, 1) and the frequencies of the 
inertial modes with [n, 1, +] = [2i, 1, +]. 

edge beams induced by the viscous oscillating Ekman layers. 
VIII. CONCLUSION 

In symmetric containers, i.e. axisymmetric containers or containers whose side walls are either parallel 
or perpendicular to the rotation axis, longitudinal libration produces an inviscid type of fluid response that 
conserves absolute vorticity. However, by friction at bottom, top and side walls, this leads to a periodic mass 
flux in boundary layers adjacent to the horizontal walls. Convergence and divergence of these Ekman fluxes 
near the edges spawn free shear layers that are inclined with respect to the rotation axis i 16 i 25 ' 33 Their angle 
is set by the ratio of the modulation frequency to the Coriolis frequency. Expulsion of Ekman mass fluxes 
takes place where the topography has the same inclination, at so-called critical slopes. In containers where 
changes in the orientation of the boundary occur suddenly (such as near the edges of the cube studied here) 
these critically sloping regions are "buried" in these horizontal edges. The free shear layers thus spawn from 
the edges, and are here termed "edge beams" . 

In a 2D approximation, these edge beams propagate along periodic paths. For paths that are short enough, 
upon a few reflections energy folds back onto itself and an inertial mode may establish itself: a vertically 
standing large-scale structure that may either be propagating or standing ( "sloshing" ) horizontally. In our 
experiment, symmetries, due to the way that the fluid is forced into motion by means of libration, allow 
the realization of only a particular subset of all possible theoretical inertial modes. Latitudinal libration 
(libration along an axis normal to the rotation axis) might provide a mechanism by means of which also the 
other, less symmetric modes can be excited, e.g. the mode [1,1,-], which consists in a single oscillating cell. 

Previous studies in cylindrical domains suggest that the resonant forcing (and, therefore, amplification) 
of a vertically standing inertial mode opens the possibility of elliptic instability (vortex breakdown). This 
should be manifest as a sudden collapse of the inertial wave mode. This has so far not been seen in our 
experiments in a cube, but is not ruled out either. Another peculiarity of the theoretical inviscid wave field 
in the cube is the appearance of many horizontal scales of motion, even for "large-scale" modes (i.e. having 
low vertical wave number and high frequency) £ These scales are required by the wave field in order to match 
the rotational (circular) to the geometrical (square) symmetries. More detailed experimental observation of 
an inertial mode is required to investigate this property. Finally, it is of interest to examine to what extent 
edge beams and inertial modes can play a role in containers whose symmetry is broken. 

Appendix A: The inviscid base flow in a librating cube 

In this section we derive the streamfunction describing the inviscid flow generated by the libration of a 
parallelepiped container of square cross section. The Euler equation in the librated frame of reference writes 

— + (« + 2fi) x u = -V(J + -Mt) x if + -|u| 2 J - -xr, (Al) 

with il(t) the total angular velocity, of vertical component given by Eq. ([T]). This equation contains three 
terms related to the (modulated) background rotation: (1) the Coriolis acceleration (last term on left-hand 
side), (2) the centrifugal force (second term on right-hand, that can be absorbed in a reduced pressure) and 
(3) the Euler force (last term on right-hand side), due to the modulation of the rotation rate^ Here u> is 
the relative vorticity vector, u> = V x u. Taking the curl of the Euler equation (|A1[) . we obtain the vorticity 
equation: 

— + u • V(w + 2H) -(<*> + 20) • Vu = -2—. 
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This equation tells that relative vorticity changes due to advection of absolute vorticity, 

u> a = ui + 217, 

(second term on left-hand side), vortex tilting by sheared motion (third term), or by changes in the rotation 
rate (right-hand side). 

Note that for forcing by libration, the driving term on the right of the vorticity equation is present in the 
vertical direction only and is independent of the vertical coordinate. Since the boundaries are not inclined 
relative to the rotation axis this motivates looking for velocity fields that are z-independent as well. In fact, 
the vanishing of the vertical velocity at the horizontal bottom and top walls suggests looking for solutions 
of zero vertical velocity. The vertical component of the vorticity uj a = uj a ■ e z (which at this point is the 
only nonvanishing component) thus satisfies 

Since the horizontal velocity field is nondivergent (d x u x + d y u y = 0), we can introduce a streamfunction 
such that (u Xl u y ) — (—d y ip,d x ip). The spatial part ^(x,y) of streamfunction, ip(x,y,t) = ^(x, y) cos(eroi), 
therefore satisfies a Poisson equation, 

Af = -2en , (A2) 

subject to the condition that in the librating frame the solid boundary at x = ±L and y = ±L acts as a 
streamline. We introduce the normalized coordinates (x, y) = (x,y)/L. The solution of Eq. (|A2[) is the sum 
of a particular solution, and a homogeneous solution, ^/j, such that A^h = 0. A particular solution is 

that matches the right hand forcing term and satisfies the boundary condition at x = ±1, but is nonvanishing 
at y = ±1. The homogeneous solution, ^h, also vanishes at J = ±1, but annihilates ^ p (x) at y = ±1. It is 
given by 

*» = -^o}j-l) dlcosHdn) , 

where 

d n = (2n+l)|. 

For convenience, we finally symmetrize the solution, by simply computing (^(x^y) + ^(y, x))/2, yielding 

fr = eQ (l 1 (.t 2 I f) 2V( 1 )» cos ( c; " ;£ ) cosh ( (i ^) +cos(d n y)cosh(d n x) \ 
°V 2 ' j-f 4cosh(d„) / 
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